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Abstract 

Organogenesis is an important process for plant regeneration by tissue or cell mass differentiation to regenerate a complete 
plant. MicroRNAs (miRNAs) play an essential role in regulating plant development by mediating target genes at 
transcriptional and post-transcriptional levels, but the diversity of miRNAs and their potential roles in organogenesis of 
Acacia crassicarpa have rarely been investigated. In this study, approximately 10 million sequence reads were obtained from 
a small RNA library, from which 189 conserved miRNAs from 57 miRNA families, and 7 novel miRNAs from 5 families, were 
identified from A. crassicarpa organogenetic tissues. Target prediction for these miRNAs yielded 237 potentially unique 
genes, of which 207 received target Gene Ontology annotations. On the basis of a bioinformatic analysis, one novel and 13 
conserved miRNAs were selected to investigate their possible roles in A. crassicarpa organogenesis by qRT-PCR. The stage- 
specific expression patterns of the miRNAs provided information on their possible regulatory functions, including shoot bud 
formation, modulated function after transfer of the culture to light, and regulatory roles during induction of organogenesis. 
This study is the first to investigate miRNAs associated with A. crassicarpa organogenesis. The results provide a foundation 
for further characterization of miRNA expression profiles and roles in the regulation of diverse physiological pathways 
during adventitious shoot organogenesis of A. crassicarpa. 
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Introduction 

With the characteristics of rapid growth, high pulp yield, high 
fiber quality and ability to thrive in degraded soils [1], Acacia 
crassicarpa has become an important economic species of the 
Fabaceae family, and has been planted for reforestation, 
reclamation of wasteland, and industrial material production on 
a large scale in Southeast Asia [2-3] . However, its recalcitrance of 
regeneration, long life cycle, and the prolonged period needed for 
evaluation of traits at maturity strongly limit traditional breeding 
programs for the species. Development of an effective propagation 
strategy is required urgently. 

Plant organogenesis in vitro is an efficient technique for large- 
scale propagation and popularization of selected genotypes. The 
technology has been extensively applied in plant genetic 
engineering, double haploid or polyploid breeding, and asexual 
reproduction of mutants or threatened species [4-5], for example. 
However, it shows strong species dependence and genotype 
specificity. For example, Arabidopsis thaliana and Mcotiana tahacum 
are easily induced to generate adventitious shoots or roots [6-7], 
whereas some plant species are recalcitrant to plant regeneration 
in vitro, such as Triticum aestivum, J^ea mays, and Gos.vypium spp [8]. 
The regeneration frequency of leguminous trees may be quite poor 
in natural habitats [9]. Therefore, a thorough understanding of the 



molecular mechanisms and gene regulatory networks involved in 
organogenesis of leguminous trees is essential in order to achieve 
improved in vitro plant regeneration and genetic transformation 
frequencies. Although important results have been obtained 
regarding the hormonal regulation of organogenesis and organ- 
related expression of genes and proteins in leguminous trees, few 
studies have focused on identification and expression of micro- 
RNAs (miRNAs) during organogenesis. 

MiRNAs are small endogenous non-coding RNA usually 
associated with gene silencing by guiding cleavage of complemen- 
tary mRNAs or suppression of translation [10]. The importance of 
miRNAs has been realized owing to their wide occurrence in all 
kinds of organisms and their important biological function in 
regulation of gene expression [11-12]. Increasing evidence shows 
that the length of 20-24 nt miRNAs plays a crucial role in 
developmental and physiological processes, including develop- 
mental timing, tissue-specific development, and stem cell mainte- 
nance and differentiation [13-14]. Plant organogenesis and 
somatic embryogenesis are the most important strategies for 
in vitro plant regeneration [15]. During the last decade, several 
reports have demonstrated the crucial roles played by miRNAs 
during somatic embryogenesis. Thus, miRNA-mediated modula- 
tion of somatic embryogenesisis is closely associated with the plant 
regeneration process [16-17]. Plant organogenesis, which is also 
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Figure 1. Morphology of different stages during plant organogenesis in Acacia crassicarpa. S1: Zygotic embryo, excised from the mature 
seeds. S2: Zygotic embryo differentiated, after two weel<s subculture. S3: Embryogenic callus, after three weeks subculture. S4: Shoot buds, after four 
weeks subculture. S5: Clusters of adventitious shoots, after 5 weeks subculture. S6: Adventitious shoot elongation. S1-S4 were cultured on MS 
medium containing 4.54 nM TDZ and 2.85 i^lVI lAA, observed using an Leica stereomicroscope. S5-S6 were cultured on MS medium supplemented 
with 2.89 |xM GA3. 

doi:1 0.1 371 /journal.pone.0093438.g001 



strongly associated with cell differentiation, undoubtedly is affected 
by miRNA repression. However, few correlative studies of miRNA 
function in the plant organogenesis process have been reported, 
especially with regard to miRNA of ^. crassicarpa. 

In order to study the regulatory role of miRNAs in organogenesis 
of ^. crassicarpa, high-throughput sequencing technology ntegrated 
with bioinformatic analysis was employed to identify conserved and 
novel miRNAs involved in adventitious shoot organogenesis of ^. 
crassicarpa, followed by prediction of ^. crassicarpa miRNA targets 
and functions. The expression patterns of one novel and 13 
conserved miRNAs during the sbc developmental stages of A. 
crassicarpa organogenesis were comparatively analyzed. This study is 
the first report focused on the modulatory roles of miRNAs during 
adventitious shoot organogenesis in A. crassicarpa. 

Materials and Methods 

Plant Materials 

The experimental materials were obtained at different devel- 
opmental stages of an adventitious bud regeneration system 
induced from mature zygotic embryos of ^. crassicarpa, as described 
by Yao et al. [18](Figure 1). Based on morphological and 
anatomical traits, six developmental stages during shoot regener- 
ation were defined and sampled. Generation of embryogenic callus 
(Figure 1, S1-S3) was induced in the dark and differentiation of 
shoot buds (Figure 1, S4-S6) occurred under light. All samples 
were frozen immediately in liquid nitrogen and stored at — 80°C. 

Extraction of Total RNA, Construction of sRNA Library and 
Solexa Sequencing 

Total RNA at the six developmental stages was extracted using 
the mirVana™ miRNA Isolation Kit (Ambion, Austin, TX, USA) 



in accordance with the manufacturer's instructions. The quality of 
the isolated total RNA was assessed by analysis with a NanoDrop 
ND-2000 spectrophotometer (NanoDrop, Wibnington, DE, USA). 
The integrity of the total RNA was monitored with a Bioanalyzer 
2100 and RNA 6000 Nano LabChip Kit (Agilent, Palo Alto, CA, 
USA). The RNA integrity number was higher than 6.0, therefore 
the isolated total RNA was suitable for library preparation. 

Equal amounts of total RNA at each developmental stage were 
pooled. About 30 |lg of total RNA was used for library 
construction and sequencing. A small RNA library was generated 
using the lUumina TruSeq Small RNA Sample Prep Kit (lUumina, 
Hayward, CA, USA) in accordance with the manufacturer's 
instructions. The purified small RNA fractions between 10 and 
40 nt were hgated with proprietary 5' and 3' adaptors separately, 
and then converted to cDNA by RT-PCR. The purified cDNA 
library was used for cluster generation on an lUumina Cluster 
Station and then sequenced with an lUumina GAIIx platform 
following the manufacturer's instructions. 

Bioinformatic Analysis of Sequencing Data 

Raw sequencing reads were obtained using Sequencing Control 
Studio version 2.8 software (lUumina) foUowing real-time sequenc- 
ing image analysis. The extracted sequencing reads were first 
processed with a proprietary pipeline script, ACGTlOl-miR v4.2 
(LC Sciences, Houston, TX, USA), for removal of 'impurity' and 
unmappable sequences resulting from sample preparation, se- 
quencing chemistry and processes, and the optical digital 
resolution of the sequencer detector. The processed raw data 
have been uploaded to Short Read Archive of NCBI with the 
accession number of SRX469994. The remaining sequences with 
lengths of 15 to 30 nt were mapped to precursor miRNAs (pre- 
miRNAs) sequences of aU plants in miRBase 19.0, for aU mature 



Table 1. 


Forward primer sequences used for qRT-PCR. 




Number 


mlRNA name 


Primer sequence (5 -3') 


1 


acr-miR159a-3p 


GCTTTGGATTGAAGGGAGCTCTAA 


2 


acr-miR162 


GCTCGATAAACCTCTGCATCCAAA 


3 


acr-miR390a-5p 


CTCAGGAGGGATAGCGCCA 


4 


acr-rriR396 


CCTTCCACAGCTTTCTTGAACTGA 


5 


acr-miR319a 


GCTTGGACTGAAGGGAGCTC 


6 


acr-miR156r 


GCGGCCTTGACAGAAGATAGAGAGCATAAA 


7 


acr-miR164a 


GGAGAAGCAGGGCACGTGCA 


8 


acr-miR166h-3p 


GCTCTCGGACCAGGCTTCATTCC 


9 


acr-miR167a 


GCGTGAAGCTGCCAGCATGATCTAAAA 


10 


acr-miR168a 


CGCTTGGTGCAGGTCGGGAA 


11 


acr-miR171q-3p 


GCGTGATTGAGCCGTGCCAATATCA 


12 


acr-miR397a-5p 


GGGTCATTGAGTGCAGCGTTGATGAAA 


13 


acr-miR398d 


GCTGTGTTCTCAGGTCGCCCC 


14 


acr-novel2* 


GTCTTTCCAATGCCTCCCATACC 


15 


5.8S rRNA 


CGCCTGCCTGGGTGTCACAA 


doi:l 0.1 371 /journal.pone.0093438.t001 
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miRNAs included in the database. The mapping sequences were 
further aligned against Acacia mangium expressed sequence tags 
(ESTs) and the genomic sequence of selected plant species, 
including Arabidopsis thaliana, Populus trichocarpa and soybean, to 
identify conserved miRNAs in A. crassicarpa. The unmapped 
sequences matching sequences in other defined databases, mainly 
comprising mRNA, RFam, and Repbase, were removed. A 
BLAST search with the remaining sequences against the selected 
genome such as Arabidopsis thaliam, Populus trichocarpa and soybean 
and ESTs of Acacia mangium was performed, and those extended 
sequences at the mapped genome positions with the propensity to 
form hairpins (as predicted with Mfold software) were selected as 
putative novel miRNAs. 

Prediction of Potential Target mRNAs 

Target prediction for the miRNAs was based on the principle of 
nearly perfect complementation between the miRNA and target 
mRNAs [19]. The identified conserved and putative novel 
miRNAs were all submitted for target gene prediction using 
TargetFinder (http:/ /hercules. tigem.it/TargetFinder.html) with 
the default parameters, using the algorithm described by Quandt 
et al [20] . The target genes were then validated using public 
resources of The Arabidopsis Information Resource (TAIR; 
http://arabidopsis.org/). Sequences with a score of less than 
4 were regarded as miRNA target genes. On the basis of their 
functions putative targets were classified using Gene Ontology 
(GO) annotations. All GO annotations were obtained from the 
TAIR database, and most (about 87%) of the targets were 
annotated. Given that many targets have more than one GO hit, 
the number of total hits can be higher than the total number of 
targets. 

Quantitative Real-time PCR 

Small RNAs of ^. crassicarpa at the six organogenesis stages were 
isolated using micro RNA Extraction Kit (BioTeKe, Beijing, 
China) following the manufacturer's instructions. With minor 
modifications, poly (A)-tailed quantitative real time PCR (qRT- 
PCR) was chosen for experimental identification. miRNA qRT- 
PCR was carried out using the NCode™ VILO'^^ miRNA 
cDNA Synthesis Kit (Invitrogen, Carlsbad, USA) and SYBR 
Premix Ex Taq^"^ (Invitrogen, Carlsbad, USA) with the small 
RNAs as the template. Investigated miRNA was amplified using 
forward primers that were designed based on mature miRNA 



sequences (Table 1) and reverse primers provided with the kit. The 
qRT-PCR amplification protocol was as follows; initial activation 
at 94°C for 2 min, followed by 40 cycles consisting of 94°C for 
20 s and 60°C for 34 s. To verify the absence of contamination, 
negative controls (no cDNA template) were carried out. The 
experiments were performed for five technique rephcates and 
three biological replicates. The data were calculated using the 
2-AA6T j^gd-jQji pi]. 5.8S rRNA cloned from A. crassicarpa was 
selected as the reference gene for normalization. 

Statistical Analysis 

The statistical analysis was performed using one-way Between- 
groups Linkage with SPSS 18.0 (SPSS Taiwan Corp., Taiwan). 

Results 

Small RNA Sequence Profile 

A small RNA library was constructed using pooled RNAs 
isolated at six developmental stages of adventitious shoot 
regeneration, and sequenced using Solexa high-throughput 
technology in a single lane. More than 10 million raw reads were 
yielded, which were filtered for adaptors and junk contaminants. 
After removal of redundant sequences we obtained 980,584 
unique sequences ranging from 1 5 to 30 nt in length, of which the 
majority was 19-26 nt long. A major peak in the abundance of 
mapped small RNAs was observed at 24 nt, which accounted for 
55.3% of the sequences (Figure 2). Those sequences were perfectiy 
mapped to pre-miRNAs of miRBase 19.0 and selected species 
genomes including Arabidopsis thaliana, Populus trichocarpa and 
soybean or Acacia mangium ESTs. The generated sequence was 
deemed a "mappable sequence". The unmappable reads were 
further searched against noncoding RNAs (rRNA, tRNA, snRNA, 
and snoRNA) deposited in the Rfam and Repbase databases and 
comprised only a small fraction (4.59%) of the counts of all 
sequences, and were removed to avoid influencing miRNA 
identification (Table 2). Finally, 2,602,396 cleaned sequences 
representing 729,967 unique reads were used for subsequent 
analysis. 

Conserved miRNA Families and Isoforms 

To identify conserved miRNAs in A. crassicarpa, we selected the 
18-25 nt mappable small RNAs for further analysis by observa- 
tion of miRNA length in miRBase. From a series of BLAST 



60% 




15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 
Length (nt) 



Figure 2. Length distribution of sequenced small RNAs. 'Redundant' represents the total number of sequences, 'unique' represents the 
number of unique sequences, nt, nucleotides. 
doi:1 0.1 371/journal.pone.0093438.g002 
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Table 2. Summary of reads from raw data to cleaned sequences of Acacia crassicarpa. 





Category 


Total reads 




Unique reads 




raw reads 


10,378,615 


1 00% 


980,584 


100% 


filtered reads 


7,315,987 


70.49% 


203,692 


20.78% 


Rfam 


449,801 


4.33% 


44,866 


4.58% 


Repeats 


26,614 


0.26% 


4,116 


0.42% 


clean reads 


2,602,396 


25.07% 


729,967 


74.44% 



Note: filtered reads: 5' adapter, 3' adapter, 5' adapter and 3' adapter joined together without insertion length with <15 nt and >32 nt and junk reads. 
doi:l 0.1 371 /journa!.pone.0093438.t002 



searches against miRBase 19.0 and Arabidopsis thaliana, Populus 
trichocarpa and soybean genomes or Acacia mangium or ESTs, we 
obtained 189 conserved unique sequences from 57 miRNA 
families allowing for sequence mismatches of not more than two 
nucleotides (Table 3). Among these sequences, the bulk of the 
miRNA families were conserved as we found orthologues of 
known miRNAs from all other plant species, such as miR156, 
miRl 66, miRl 67 and miR396, whereas some miRNA families were 
less conserved and were present only in some of the plant species, 
such as miRl 70 (conserved in Arabidopsis only), miR6478 and 
miRl 448 (weakly conserved in Populus trichocarpa). Eight miRNA 
families were species-specific to the leguminous model plant 
soybean or Medicago truncatula, comprising miR2118, miR4415, 
miR5368, miR4995, miR5281, miR2592, miR5740, and miR5256. 

The majority of conserved miRNA families contained members 
that differed in no more than two nucleotides, which are regarded 
as miRNA isoforms [22] . Interestingly, the number of isoforms in 
differently conserved miRNA families varied considerably. Four 
families [miRl 56, miRl 66, miRl 67, and miR396) contained the 
highest number of multiple isoforms with 13, 11, 10, and 10 
members, respectively. Eighteen miRNA families (e.g., miR165, 
miR170, miR393, miR399, and miRlSll) contained only one 
member. 

The abundance of known miRNA families was also monitored. 
Interestingly, a notable divergence in expression frequency among 
miRNAs during shoot regeneration was observed. For example, 
miRl 59, miRl 66, and miR396 were detected most frequently 
with the counts of 15,425, 11,560, and 15,290, respectively; 
miR156, miR157, and miR319 showed moderate abundance with 
counts of 4355, 3566, and 4020, respectively; and miRNAs such 
as miR2592, miR5256, miR5740, and miR6485 showed low 
abundance (Table 3). 



Putative Novel miRNAs 

The sequences unmapped to the pre-miRNAs in miRBase 
19.0 were chosen to identify putative novel miRNAs. The 
following criteria were applied: (1) the unique sequences were 
mapped to the selected species genomes; (2) the extended sequence 
at the mapped genome location had the propensity to form a 
hairpin. Taken together, a total of 7 novel miRNA candidates 
from 5 families were obtained with a minimal folding free energy 
index (MFEI) for the hairpin structure of the miRNA precursor of 
less than 0.7 (Table 4). These miRNA candidates commonly 
displayed a strong bias toward a 5' U at the first nucleotide 
position (Figure 3) and stem-loop structure precursors (Figure 4), 
which was consistent with conserved miRNAs [23]. Some previous 
studies reported that a standard for discrimination between high- 
confidence microRNAs and fragments of other RNAs in deep 
sequencing data was the detection of miRNA* sequences [23]. 
Among the miRNA candidates in the present study, four novel 
miRNAs were identified with their complementary miRNA* 
belonging to two miRNA families, which implied that these 
putative miRNAs were most likely novel to A. crassicarpa, and 
several potentially novel miRNAs might be specific to A. crassicarpa 
or Fabaceae species. Of these potentially novel miRNAs, the copy 
number of several mLRNA*s was low. This feature is consistent 
with miRNA* degradation during the miRNA generation 
procedure [24], so counts of miRNA*s were not included in the 
most recent criteria for annotation of plant miRNAs. In the 
present study, acr-novel2* showed higher counts than the corre- 
sponding novel miRNAs. This finding was in accordance with the 
recent discovery by Zhang d al. [25] that miRNA* sequences are 
abundantly expressed. 
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Figure 3. First nucleotide bias of identified miRNAs in Acacia crassicarpa. (A) Conserved nniRNAs, (B) novel miRNA candidates. 
doi:1 0.1 371/journal.pone.0093438.g003 
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Table 3. Identification of conserved miRNA families in Acacia crassicarpa with corresponding isoforms. 



miRNA family 


Gene name 


Length 


Sequence 


Counts 


Reference miRNA 


MIR1 '^fi 

IVIln 1 JO 


3cr~miR1 56d 


20 




1 


11 iLi 11 im 1 DDd 




3cr-miR1 56b 


20 


1 if^ArAi^AAi^Ar^Ai^i if^Af^rai i 


g 


rr^t miR 1 '^AK 
CCa min 1 DDD 




3cr-miR1 56c 


23 


1 ir;A("A(^AAi^A(^Ar;i \r:.^(^r^rrrh 


1 
1 


ghr-mlRI 56c 




3cr-miR1 56d 


20 


1 if^ArA(^AAi^A(^Ar;i if^Af^("Ar~ 


Do'l- 


Qm3-miR1 563 




3cr-miR1 56g 


21 


1 II i(^Ar'Ar;AAf^Ar;Ar;("(^A(^r'Ar' 

U Uu AL-AuHAuAuHu^uAuL-AL 


an 


5DI min 1 DDc 




3cr-miR1 566-p3 


21 


r-.r\ \r\ in iai ii in ii in if^i ii^ai \r 


44 


LLC Ivlln 1 DOe pD 




arr mlDI ^Af 

acr mi n i dot 


21 


1 if^ArAi^AAi^Ai^Ar;! 1 A f^(" A rr" 

UvjAV-AVjAAuAuAuUoAvjV-AH— 




n la mi n i DDa 




acT'iniRI 56q 


1 8 


ArAf^AA(^A(^A(^l ]CAcrAr 

A*— AOAAVJ AUA\JUOA^J'_A^_ 


5 


ath-miRI An 
d Lr 1 [ T II n 1 DOy 




acr~miR1 56h"5p 


21 


l 11 if^ArAf^A Af^ArAr;Ar;A(^rAr 

U UUAV-AVJ/AAVJA\^AU AUAVJV_A\^ 


1 


diy 1 1 iin 1 DDI 1 Dp 




acT'iniRI 56i"3p 


22 


r;ri irAn ii in in ii ii in \ci irAi ir 


5 


mtr-WIRI '^Ai-i-i^ 

1 1 III IVll n 1 DDI p J 




acr-miRI 56j 


21 


1 II if^Ar^chhc^f^f^f^i ir;Ar;r"A(" 

U UuAV-AuAAuAuuuUuAuLA^ 




mtf-mlRI 56j 




acr-miRI 56j-3p 


22 


i^n irAn ii in in ii ii in ii^i itai ii i 

Ijl^UV-A'^UULUl^UUUl.wUuU'—AUU 


2 


/-mo MIRI'^Ai 
Cme min 1 DDJ pD 




3/-r m!R1 i^AL- 
dCr min 1 DDK 


21 


1 11 If^ATAf^A Al^Af^ Ar;A(^ A(^r~Af~ 
UUuAL-AuAAuAuAvjAuAuLAL 


1 7 


yma min i ddd 




-if-r m!R1 a! 
aCr min 1 DDI 


21 


iiiir;ArAr;AAf^Ai ii^r;Ar;Ar;f"Ar" 

U UuAV-Au AAvjAUIju AoAuL AL 




nti- miR 1 ^Al 
piC min 1 DDI 




acr-miRI 56m 


21 


iiiif^ArAr^AAf^Ai iAr;A(^ Af^rAr" 

UUuA^-AlJAAv□AUAv□AuAuLAL 


3530 


gm3-miRl 56c 




acr-miRI 56n 


21 


1 II i(^ ATAi^ A Af^Ai 1 A A A r" 

U UuAV-AuAAoAUAoAuouLAL. 


■J 


mtr-miRI 56q 




3cr-miR1 56o-3p 


22 


cri irAn lAn in ii ii in it^i irAf^i i 

U V- AV- U Al^ U \^ U U U 1^ U U U Vi. A\J U 


1 


nma-MIRI "lAn-n^ 
y 1 lid ivim 1 Dou pD 




acr-miRI 56p 


21 


1 11 If^ArAl^A Af^AI lAf^nt^Af^rAr 

U U\J A>_AVJ AAOAU AOUOAVJl^ Al^ 


4 


Qm3-miR1 56p 




3cr-miR1 56r 


21 


1 11 l(^ArAr;A A(^AI IA(^A(^ Af^rAI 1 

U UVJ AV-AVJ AAVJAU AVJAU AU*— AU 


58 


Qm3-miR1 56r 




acr-miRI 57d-3p 


21 


i^n in in iai if^n ii in if^i I("ai \r 


g 


ath MIR1 ^~lA 
ain Ivlln 1 D/u pD 




3cr-miR1 57d-5p 


20 


1 ir;A("A(^AAi^Ai iAr;Ar;Af^("Ar~ 

UuAV_AuAAuAUAuAuAuV_AL 




ain min i d /a 


(vUn 1 Dy 


acr-miRI 59a-5p 


1 9 


r;Af^n irn ii ir^AAi^i ittaai i 

oAoV-UH-U UoAAuUl-l-AAU 


2 


Qm3-miR1 593-5p 




3cr-miR1 59a-3p 


21 


1 II 11 ir;i^Ai II i(^AAi^(^(^A(^n in ia 

U UUuuAU UuAAljvjuAVj>.wULUA 


1 5030 


gm3-miR1 593-3p 




acr mi n i jyo 


21 


r;Af^n irn ii ii^AAr;i irrAAi ia(^ 

oAoV-ULLU UuAAuUl-l-AAUAu 


1 o 


pic IVll n 1 DyD pD 




acr min i Dye 


1 g 


1 II II ir;r;Ai ii ir^AAi^r^i^Ai^n ir 

UUUuuAU UuAAuuuAuLUL 




ain min 1 Dye 




dv-i 1 Mini Dyu 


1 8 


1 II If^f^AI 11 If^AAf^f^f^Af^n \C 

UVJVJUAU UVJAAVJVJVJAVJl^U*— 


g 


pLi. mm 1 Dyu 




3cr-miR1 59g 


20 


1 11 II If^f^AI 11 \CA!I^CCCACri in 1 

U UUUUAU UUAAUVJVJAVJ*— Ul^U 


195 


dLii 1 Mini Dyu 


JVlIn 1 DU 


3cr-miR1 60a-5p 


21 


1 ir;m ii^i^n irrn ii^i iai ii^rrA 


1 

1 J J 


gm3-miR1 603-5p 




3cr-miR1 603-3p 


21 


i^rr;i IAI if^Ai^f^Ar^r"f"AAr;rAi ia 

U^uUAUuAuuAuLLAAuLAUA 


-J 


gm3-miRl 603-3p 




acr min i dud 


21 


rci ir;r;Ai ii^r;r(^i iai ir;Af^i^Ar;r 

HjULjuMUuuL-uUAUuAuuAuV- 


1 


ath MIRIAHa 
dTM Ivlln 1 OUa pD 


Ivlln 1 OZ 


3cr-miR1 62 


21 


1 irr;Ai iAAAr~n in ii^tai ir"r~AA 

UV-UAUAAALLUV-UuLAULLAA 




3311 miD1 AD 

aau min i dz 




acr Ivlln 1 OZ Dp 


21 


f^r;A(^(^r"Af^rf^(^i ii itai irr^Ai \r 

u o Au u Ao l_ ou U U l_ A U l_ o A U L 


D 


aai 1 ^^IR1 AT nC; 
aau Ivlln 1 dz pD 


JVI 1 n 1 D'f 


acr-miRI 64a 


21 


1 ii^(^A(^AAi^rA(^r;i^rAri^i ir;r"A 

UvjuAuAAux^Auuljl-AV-oUuLA 


AO 


gm3-miR1 643 




3cr-miR1 64c-5p 


21 


1 if^r^Ai^AAi^rAi^i^f^rArr^i \crr 

UouAuAAul-AuVjljl-A'—oUuLL 




al\/ miRIAAr CTn 

aiy min i o^c op 


MIRIAM 
ivim 1 oj 


3cr-miR1 65a-3p 


21 


1 \CC-,CiACC [ii.CCC\ II irAi \ccccc 

U \J U A 1^ 1^ A U U U U A U 1^ 1^ 1^ l^V- 


7 


diy 1 Mini DDd Dp 




3cr-miR1 65a-5p 


23 


rif^AAi ir;i ii ini in \cc-,c\ irriAt^nAi i 


1 


aK/-miR1 A^;4-'^ri 
diy 1 Mini DDd Dp 


^/l 1 R 1 A A 
J VI 1 n 1 DO 


3cr-miR1 66a 


21 


1 \rr-tiArr tiCr-r\ w ir~Ai ii \rrrr 

U\^uAAV_^Auu^UULAUU\^^\^^ 


g 


rco-miRI 663 




3cr-miR1 66a-3p 


21 


1 irf^(^Af"f"Ar;i^n ii ir~Ai ii \rrrr 

U \_ vju ALL Au U U A U U 


9042 


Qm3-miR1 663-3p 




3cr-miR1 663-5p 


21 


i^i^AAi ii^i 11 if^i in ii^i^n \rr-Lr-c 

uuAAUuUUuULUuuHJLuAuLl 


1 37 


gm3-miR1 663-5p 




acr-miRI 66b-5p 


21 


i^f^AAi ii^i II \r-,\ in ii^r;i ii \rr-Ar-c 

uvjAAUuUUuULUuvjUUHjAuu 


1 


os3-miR1 66b-5p 




3cr-miR1 66c-5p 


20 


i^r;AAi ii^i \rr-\ in ii^r;i ii irr^Ai i 


1 


i-\ti- ^/llR1AAcl 
pic Ivlln 1 DOc pD 




3cr-miR1 66d-3p 


21 


1 \rr-r-,crr tiCcr\ w itai ii \rrrr 


D 


mtr mlDIAA/^ 

mil min 1 Dou 




a/^T-mlD 1 AAH-^ r\ 
ax-t I Mini Ouu Dp 


21 


(^(^AAi \Ci\ II if^rn \c-iC-,c\ \ccacc 


1 


Mill IVIin 1 DDU pD 




3cr-miR166e-5p 


21 


GGAAUGUUGGCUGGCUCGAGG 


16 


gm3-MIR166e-p5 




3cr-miR166f-3p 


20 


UCUCGGACCAGGCUUCAUCC 


1 


bna-MIR166f-p3 




3cr-miR166g-5p 


21 


GGAAUGCUGUCUGGUUCGAGA 


1 


ptc-MIR166g-p5 




3cr-miR166h-3p 


21 


UCUCGGACCAGGCUUCAUUCC 


2271 


gm3-miR166h-3p 




3cr-miR166i 


23 


UCGGACCAGGCUUCAUUCCCCUU 


3 


cme-miR166i 




3cr-miR166i-5p 


21 


GGAAUGUCGUCUGGUUCGAGA 


n 


gm3-miR166i-5p 




3cr-miR166j 


21 


UCUCGGACCAGGCUCCAUUCC 


1 


ptc-miR166p 
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: , 

3cr~miR1 66l"3p 


21 


U \J r\ l^/A VJ U U U r\ U U U V_ V_ V_ 


10 


\j \ \ Id IVI in 1 ooi p J 




3cr-miR1 66rn 


1 9 


(^r:.Arr A(^r-.r\ ii ir~ai ii \rrrr 


ZD 


Qma-miRI 66m 




3cr-miR1 66p 


1 g 


1 \rr-r-,iirr tLr-,cr\ ii ir~Ai ii \rr 




gma-miRI 66p 


1 VI 1 n ID/ 


3cr-miR1 67a 


21 


1 \r-LLcr\ \crr tiCr ii\ \r-L\ in ia 


OH 


Qma-miRI 67a 




acr-miRI 67a-3p 


23 


r-r-\ ir~ai ir^r'i ir;i ir;ArA(^m i("Ar~i i 


c 
J 


nma ^/llR1A7a r-\5 

y ma IVIln i o/a pj 




acr-miRI 67b-3p 


21 




f. 
D 


CSI IVI! n 1 O/ D pj 




acr-miRI 67b-5p 


23 




•3 

J 


mir min 1 D/DOp 




acr-miRI 67c 


21 


1 ir;AA(^ri if^rrAf^rAi ir;Ai in ir 




vvi-miRI Q7c. 




dUI 1 1 lln 1 D/ U 


22 




266 


Qma-miRl Q7q. 




acr-miRI 67g 


22 


1 if^AArin if^rrAnrAi ir^Ai in ii ia 




Qma-miRI 67g 




acr-miRI 67G-3p 


20 


r;Ai i("Ai \r-,\ \r-r-,r iir-\ w ii \rtirr 


4 


nma ^/llR1fi7o 

y ma IVIln 1 D/ e pj 




acr min i d/t 


23 


1 if^AAr;n \crr t^cr t^x if^Ai in ir^Ar" 

UvjAAuLUuI^I^AuI.wAUvjAU'—UuHL 




dpr-miRI 67c 




acr-miRI 67g 


22 


1 ii^AAf^n ir^rrAf^rAi ir;Ai in ir^r^ 

UuAAuLUuV-V-HuLAUuAU^UIjLi 


1 n 


gma-mlRI 67c 


MIQ1 
iVIln 1 Dd 


acr min i ooa 


21 


1 irr;n ii \r-,r-\ \r-,r iir-r-\ \rr-r-,ctLti 




yma-miRI 68a 




acr min i ooa op 


21 


rrrr-rr\ ii if^rAi ifAAn ii^aai i 

LLLuV-^U UuLHUV-AAHJuAAU 




nms A^IRIARs rv^ 

yma IVIln i ooa pj 




acr min i dod 


21 


1 II \r-,r\ II ir;f^i \r-,r tLCCx xrcccti l 

U Uu^wUUIjoUuLAouUHjuuAA 


2 


mir min i Doa 




a\-i I Mini uou 3p 


22 


cccccc\ 11 ir;rAi irAAn if^AAi ii i 


g 


1 iLd iviin 1 Dod |J.3 


iviin 1 oi? 


ciLr 1 1 lin 1 D^d 


21 


rAr;rrA Af^i^Ai iriAn ii \cccc-,c-, 


5 


yma-miRI 69a 




acr-miRI 69a-5p 


21 


A! if^rAnrrAAf^f^Ai ir^An ii \c-,c 


1 


pLL ivii n 1 oyd pj 




acr-miRI 69a-3p 


23 


ccr Hex II if^i \rr\ w if^i^n iAr~Af"i i 


1 


ath MIRIAQa 

ain IVIln 1 oya p J 




acr-miRI 69aa 


21 


r"Ar;("("AA(^AAi ir;An ii ir^i \rr-c 

LAvjV-V-AAuAAUuALUUuU^UU 


1 

I 


pic min 1 oyaa 




acr mi n i oyu op 


20 


Ar;rrAa(^(^Ai ir^An ii if^rrr^i^ 

Aljl_V_AAuVjAUoA>.wUUo'_l_ou 


2 


aiy m In 1 oyaop 




acr mi n i oyo op 


21 


ccr L iir-,\ II if^i in ii ii ii^f^n iai ir^ 


■3 


yma IVIln i oya pj 




acr-miRI 69g 


24 


1 if^Af^rrAAf^f^Ai if^An ii \r-,rrccrr 

UvjAuV-V-AAooAUoAI-U UVjl_l_uljl_l_ 


C 

J 


yma-miRI 69g 




acr-miRI 69j-5p 


21 


Ai^rrAAf^AAi ir;An ii \r-rrr-r-.t^ 

rtUV-^AAuAAUuAl-UUu^LuuA 


4 


/-ima miDIAQi 

yma min i oyjop 




cit^l f Mini D^ls 


21 


TAf^rrA Af^AAl If^An )l \CCCC-,C-, 


49 


nma-miRI f,Q]e 
y 1 lid 1 Mini □? i\ 




acr-miRI 69m 


24 


1 K^Af^rrAAf^f^Ai if^An ii if^rrrif^rA 

UVJAVJ\^V_/AAU\JAUUAV-.U UU*— V_UU\-.A 


16 


vvi 1 Mini oy 1 1 1 




acr-miRI 69 p 


1 9 


1 i(^A(^rrAAr;(^Ai jf^Ari ii \gc 

UOAUV-V-AAVJOAL/UAV-U UU*— 


1 


yma-miRI 69p 




acr min i oyr 


22 


1 ii^A(^("("AAr;AAi i(^An ii \r-rrcc 

Uv3AuL.i_AAuAAUljAHJUuL.V_u\J 


ifi 


cmG-miRI 69r 




acr min i oyz 


22 


rAf^rTA Af^AAi i(^An II \r-r r c-cr 

LAui-.L.AHuAAUuAHJUlj\^l^v3UL 


2 


pic min 1 oyz 


^yl 1 R 1 7 1 

IVIln 1 / 1 


acr min i / i a 


21 


1 II \r-,i^r-rrr-{ if^rrAAi iai irArr^ 

UUuAuLUjUuLLAAUAU^A^U 


■J 


yma-miRI 71a 




acr-miRI 71 c-5p 


21 


Af^Ai IAI II ir;f^i ii^f"r;f^i II irAAi ia 

AuAUAU UuoUuHjoULH-AAUA 


1 


al\/ miR171r CTn 

aiy min 1 / 1 cop 




acr-miRI 71 b-3p 


21 


1 ii^A(^rrr;AAi itaai iai itati \r 

UvjAuV-V-UAAUV-AAUAULAl-Ul— 


Q 


yma-miRI 71 b-3p 




acr-miRI 71 G-5p 


22 


1 iai II ii^i^rn \r-,c\ ii irAn itai^aa 

UAUUuVj*-LUuVjUU'_AI_U'_A>j]AM 




mtr ft/IIRI 71 nC: 

mir IVII n 1 / 1 c p J 




acr-miRI 71 d-3p 


22 


1 11 \G[^CiCcc\ if^rrAAi iai irArAA 

U UUAU*— V^UUU*— V^AAUAU\-.AV_AA 


1 


zma-miRI 71 b-3p 




aLr 1 1 II n 1 / 1 lOp 


21 


AriAl iai 11 irit^i lArrif^i ii irAAi ir 

AUAUAU UUVJUAl— VJVJUUl— AAU<— 


5 


y 1 1 Id IVI in 1 / 1 1 p J 




al_l [Mini / IJOp 


21 


1 IAI 11 if^f^m \CiC-,\ 11 irAn irAf^A 


12 


Qma-miRI 71j-5p 




acr-miRI 71 k-3p 


21 


1 11 i(^A(^("rr;rf^rrAAi iai ir"An i 

U UuAu>.wHj^vjV_V_AAU AULA\-.U 


o 
0 


gma-miRI 71 k-3p 




acr-miRI 71 m 


21 


1 11 \ct^r-rrr-rr-\ itaai iai in ir~A 

UULiAul_^_uL.uUl_AAUAU\-.Ul-A 


1 
1 


rtm^a mlD171m 

yma min i / i m 




acr-miRI 71 c|-3p 


21 


1 \r-L\ II \c tLcr r r-,\ \crrLti\ iai ir" 

UvjAUUuAvjV-V-VjUuI-V-AAUAUL 


1 1 


nma ^/llR171n 

yma IVIln i / i q p5 


MIR 1 71 
iVI 1 n 1 / Z 


acr-miRI 72a-3p 


21 


Ai^AAi ir'i 11 ii^Ai ir;Ai ir;n ir^rAi i 

AuAAULU UuAUuAUuV-UuLAIJ 


J 


al\/ miR17'5a '^r\ 

aiy min i / zaop 




acr-miRI 72c 


21 


GGAAUCUUGAUGAUGCUGCAC 


49 


Qma-miRI 72c 




acr-miRI 72d-5p 


21 


GCGGCAUCAUUAAGAUUCACA 


7 


ptc-iVllR172d-p5 




acr-miRI 72c-3p 


21 


AGAAUCUUGAUGAUGCUGCAG 


45 


aly-miR172c-3p 




acr-miRI 72g 


22 


AGAAUCUUGAUGAUGCUGCAGU 


7 


cme-miR172e 


MIR2111 


acr-miR21 11a 


21 


UAAUCUGCAUCCUGAGGUUUU 


1 


ptc-miR21 11a 




acr-miR21 lla-Sp 


21 


UAAUCUGCAUCCUGAGGUAUA 


8 


aly-miR21 1 la-5p 




acr-miR2111b 


21 


AAUCUGCAUCCUGAGGUAUAG 


26 


gma-miR21 1 1 b 


MIR2118 


acr-miR2118a-3p 


22 


UUUCCGAUUCCACCCAUUCCUA 


912 


gma-miR2118a-3p 
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MIR2590 


acr-miR2590d-5p 


18 


UCAGCGUGGUCGGAAAUC 


1 


mtr-MIR2590d-p5 


MIR2592 


acr-miR2592bj-5p 


1 8 


AUUCCCACUGUCCCUGUC 


3 


mtr-MIR2592bj-p5 


MIR2631 


acr-miR2631 -3p 


1 9 


UUCAAUUGUAAAAUUUUGU 


2 


mtr-MIR2631 -p3 


MIR2655 


3cr"miR2655n"5p 


21 


UUUGAUCCUUUCUGUAAAUUU 


2 


mtr-MIR2655n-p5 




acr"miR26550"3p 


21 


UAAAUUUACAGAAAGGAUCAA 


2 


mtr-MIR2655o-p3 


MIR291 1 


acr-miR291 1 


23 


GCCGGCCGGGGGACGGACUGGGA 


314 


han-miR291 1 




acr-miR291 1 -5p 


23 


CCCGAACCCGUCGGCUGUCGGCG 


500 


han-MIR291 1 -p5 


MIR2916 


acr-miR2916-5p 


24 


CUAGUCUCAACCAUAAACGAUGCC 


42 


peu-MIR2916-p5 




acr-miR2916-3p 


18 


GCGGAUGUUGCUUUUAGG 


31 


peu-MIR2916-p3 


MIR319 


acr-miR31 9 


25 


UUUGGACUGAAGGGAGCUCCUAAUU 


2 


tcc-miR31 9 




acr-miR31 9a 


20 


UUGGACUGAAGGGAGCUCCC 


3259 


grna"rniR3 1 9a 




acr-miR319f 


1 9 


UUGGACUGAAGGGAGCUCC 


43 


gma~miR3 1 9g 




acr-miR319g 


21 


UUUGGACUGAAGGGAGCUCCLI 


689 


gma-miR3 1 9g 




acr-miR319h-5p 


21 


AGCUGCUGACUCAUUCAUUCA 


12 


gma-MIR31 9h-p5 




acr-miR31 9i 


21 


U UUGGACUGAAGGGAGCUCCC 


1 7 


ptc-miR31 9i 




acr-miR319n 


21 


UUUGGACCGAAGGGAGCUCCU 


1 


gma"rniR3 1 9n 


MIR390 


acr-miR390a 


22 


AAGCUCAGGAGGGAUAGCGCCU 


3 


crnG-miR390a 




acr-miR390a-3p 


21 


CGCUAUCCAUCCUGAGUUUCA 


1 


gma-rniR390a-3p 




acr-miR390a-5p 


21 


AAGCUCAGGAGGGAUAGCGCC 


398 


gma-miR390a-5p 




acr-miR390g-3p 


21 


CGCUAUCUAUCCUGAGUUUUA 


3 


gma-MIR390g-p3 


MIR393 


acr-miR393a 


22 


UCCAAAGGGAUCGCAUUGAUCU 


9 


gma"rniR393a 




acr-miR393b-3p 


21 


UAGGAUCAUGCUAUCCCUUUG 


1 


gma-MIR393b-p3 




acr-miR393i-3p 


21 


AUCAUGCGAUCUCUUAGGAAU 


2 


gma-MIR393i-p3 


MIR394 


acr-miR394a-5p 


20 


UUGGCAUUCUGUCCACCUCC 


69 


gma"miR394a"5p 


MIR395 


acr-miR395a 


21 


UGAAGUGUUUGGGGGAACUCC 


104 


gma"rniR395a 




acr-miR395a-5p 


18 


GUUCCCCUGAACACUUCA 


2 


gma-MIR395a-p5 




acr-miR395b 


21 


CUGAAGUGUUUGGAGGAACUC 


2 


cca-miR395b 




acr-miR395c 


21 


CUGAAGUGUUUGGGGGAACUC 


24 


ptc-miR395b 




acr-miR395d 


21 


CUGAAGUGUUUGGGGGAACUU 


1 


gma-miR395d 




acr-miR395g-5p 


21 


AGUUCCUCUGAACACUUCACC 


63 


gma-MIR395g-p5 




acr-miR395h-5p 


21 


AGUUCCUCUGAACGCUUCAUU 


1 


gma-MIR395h-p5 




acr-miR395i 


20 


UGAAGUGUUUGGGGGAACUC 


3 


gma~rniR395i 


MIR396 


acr-miR396 


22 


UUUCCACAGCUUUCUUGAACUG 


7 


dpr-miR396 




acr-miR396c 


21 


UUCCACAGCUUUCUUGAACUG 


1 3334 


aaU"miR396 




acr-miR396-3p 


21 


GUUCAAUAAAGCUGUGGGAAG 


167 


aau-MIR396-p3 




acr-miR396a-5p 


22 


AUUCCACAGCUUUCUUGAACUG 


56 


gma-rniR396a-5p 




acr-miR396a-3p 


21 


AUUCAAGAAAGCUGUGGAAAA 


13 


csi-MIR396a-p3 




acr-miR396b-5p 


21 


UUCCACAGCUUUCUUGAACUU 


261 


gma-miR396b-5p 




acr-miR396b-3p 


21 


GCUCAAGAAAGCUGUGGGAUA 


29 


gma-miR396b-3p 




acr-miR396e 


22 


UUCCACGGCUUUCUUGAACUGC 


32 


cm2"miR396G 




acr-miR396f 


21 


UUCCACGGCUUUCUUGAACUG 


1 388 


ptc-miR396f 




acr-miR396g-5p 


19 


UUCCACGGCUUUCUUGAAC 


2 


ptc-miR396g-5p 




acr-miR396j-5p 


21 


UUCCACAGUUUUCUUGAACUG 


3 


gma-MIR396j-p5 


MIR397 


acr-miR397a 


21 


AUUGAGUGCAGCGUUGAUGAA 


28 


gma-miR397a 




acr-miR397a-5p 


21 


UCAUUGAGUGCAGCGUUGAUG 


12 


aly-miR397a-5p 




acr-miR397b-3p 


21 


AUCGACGCUGCACUCAAUCAU 


43 


gma-miR397b-3p 




acr-miR397c 


21 


CAUUGAGUGCAGCUUUGAUGA 


1 


ptc-miR397c 


MIR398 


acr-miR398 


22 


UGUGUUCUCAGGUCGCCCCUGU 


3 


nta-miR398 




acr-miR398a 


23 


UGUGUUCUCAGGUCGCCCCUGUU 


1 


bdi-miR398a 



PLOS ONE I www.plosone.org 8 April 2014 | Volume 9 | Issue 4 | e93438 



miRNAs during Organogenesis of Acacia crassicarpa 



Table 3. Cont. 



miRNA family 


Gene name 


Length 


Sequence 


Counts 


Reference miRNA 




acr-miR398d 


21 


UGUGUUCUCAGGUCGCCCCUU 


4 


gma-miR398a 




acr-miR398b 


21 


UGUGUUCUCAGGUCCCCCCUG 


1 


ath-miR398b 




acr-miR398c 


21 


UGUGUUCUCAGGUCGCCCCUG 


181 


gma-miR398c 


MIR399 


acr-miR399a 


20 


UGCCAAAGGAGAGUUGCCCU 


1 


gma-miR399a 




acr-miR399b 


21 


UGCCAAAGGAG AAU UGCCCUG 


2 


vun-miR399b 




acr-miR399d 


21 


UGCCAAAGGAGAGUUGCCCUU 


8 


cme-miR399d 


MIR403 


acr-miR403-3p 


21 


UUAGAUUCACGCACAAACUCG 


975 


aly-miR403-3p 




acr-miR403a 


19 


UUAGAUUCACGCACAAACU 


3 


gma-miR403a 


MIR408 


acr-miR408-5p 


21 


CUGGGAACAGGCAGAGCAUGG 


44 


ptc-miR408-5p 




acr-miR408a-3p 


21 


AUGCACUGCCUCUUCCCUGGC 


157 


gma-mlR408a-3p 




acr-miR408d-5p 


21 


GCUGGGAACAGGCAGAGCAUG 


19 


gma-MIR408d-p5 


MIR4414 


acr-miR4414b 


21 


UGUGAAUGAUGCAGGAGCUAA 


5 


mtr-miR4414b 


MIR4415 


acr-miR4415a-3p 


21 


UUGAUUCUCAUCACAACAUGC 


305 


gma-miR441 5a-3p 


MIR4416 


acr-miR4416c-3p 


21 


AUACGGGUCGCUCUCACCUGG 


2 


gma-MIR4416c-p3 


MIR482 


acr-miR482-5p 


21 


GGAAUGGGCUGUUUGGGAAGC 


1 


pvu-miR482-5p 




acr-miR482a-5p 


20 


GGAAUGGGCUGUUUGGGAAG 


7 


mdm-miR482a-5p 


MIR4995 


acr-miR4995-5p 


20 


UGGCUUGGUUAAGGGAACCC 


5 


gma-MIR4995-p5 


MIR5015 


acr-miR5015-5p 


24 


UUUUGUUGUUGUUGUUAUUAUGUU 


1 


ath-MIR5015-p5 


MIR5030 


acr-miR5030-3p 


20 


UUCCUGAAGAACAAAUUGUU 


1 


gma-MIR5030-p3 


MIR5139 


acr-miR5139 


18 


AACCUGGCUCUGAUACCA 


282 


rgl-miR5139 


MIR5248 


acr-miR5248 


21 


UUUUUAGUUGGCAUGCAUUUA 


2 


mtr-miR5248 


MIR5256 


acr-miR5256-5p 


21 


UUGUAAGAUUAAAAUGGUUGA 


3 


mtr-MIR5256-p5 


MIR5257 


acr-miR5257-5p 


21 


UCAUCAAGAACAAGUAGAACU 


1 


mtr-MIR5257-p5 


MIR5261 


acr-miR5261 


21 


UUCAUUGUAAAUGGCUUUGGC 


1 


mtr-miR5261 


MIR5281 


acr-miR5281 b 


24 


CUUAUAAUUAGGACCGGAGGGAGU 


6 


mtr-miR5281b 




acr-miR5281c-5p 


23 


UAU AU AU UACU ACCUUCGGUCCU 


1 


mtr-iVllR5281c-p5 


MIR5282 


acr-miR5282-3p 


24 


GCAAAAUUUUGUGACGGAAUUAGU 


1 


mtr-iVllR5282-p3 


MIR530 


acr-miR530a 


22 


UGCAUUUGCACCUGCACUUUAA 


1 


gma-miR530a 


MIR5368 


acr-miR5368-3p 


19 


GAUACCACUCUGGAAGAGC 


96 


gma-MIR5368-p3 




acr-miR5368-5p 


20 


CCUGGGAUUGGCUUUGGGCC 


16 


gma-MIR5368-p5 


MIR5740 


acr-miR5740-5p 


21 


UUGAUUUGUAUCUGUUUGGAU 


3 


mtr-iVllR5740-p5 


MIR6478 


acr-miR6478 


25 


CCGACCUUAGCUCAGUUGGUAGAGC 


190 


ptc-miR6478 


MIR6485 


acr-miR6485-5p 


23 


UCGGCAGAUUUGGAUUCCUAUAU 


3 


hbr-MIR6485-p5 


MIR827 


acr-miR827 


18 


UUAGAUGGCCAUCAACGA 


1 


pUc-miR827 


MIR828 


acr-miR828 


22 


ACUUGCUCAAAUGAGUAUUCCA 


1 


gma-miR828a 


MIR858 


acr-miR858a 


21 


UUUCGUUGUCUGUUCGACCUU 


12 


ath-miR858a 




acr-miR858b 


21 


UUCGUUGUCUGUUCGACCUUG 


7 


ath-miR858b 


MIR869 


acr-miR869-5p 


18 


CAACUCCAGGAUUGAACC 


1 


ath-MIR869-p5 


MIR1310 


acr-miR1310-3p 


20 


AAACUUUAAAUAGGUAGGAC 


66 


han-MIR1310-p3 


MIR1448 


acr-miR1448-3p 


20 


UUUCCAAUGCCUCCCAUACC 


14 


ptc-MIR1448-p3 


MIR1511 


acr-miR1511 


18 


AACCAGGCUCUGAUACCA 


11 


gma-miRI 51 1 



doi:l 0.1 371 /journal.pone.0093438.t003 



Putative Target Prediction for Conserved and Novel 
miRNAs 

Putative targets of identified conserved and novel miRNAs were 
predicted by BLAST using niRNA sequences from A. thaliana [22] 
using TargetFinder. On the basis of completely complementarity 
between the miRNA and target mRNA [26], we screened 237 
genes for feasible targets with a total score not higher than 3 
(Table 5). Among the targets predicted, the most frequent were 



transcription factors, including a SBP transcription factor, MYB 
transcription factor, NAG transcription factor and GRF transcrip- 
tion factor, and other putative target genes involved in a range of 
physiological or metabolic processes, for example functional 
proteins such as protein kinases, laccase, and protein phosphatase. 
Many reports have demonstrated that target genes of conserved 
miRNAs are consistent among the majority of plant species [27], 
and this was the case in the present study. Putative target genes of 
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Figure 4. Stem-loop structure of partial novel miRNA precursors. (A) acr-novell and acr-novell*, (B) acr-novel2 and acr-novel2*, (C) acr- 

novel3, (D) acr-novel4,(E) acr-novel5. Novel miRNA sequences were noted by red and yellow shadow. 

doi:10.1371/journal.pone.0093438.g004 



miR156/ 157 mainly encoded a SBP transcription factor family 
protein and SQUAMOSA promoter binding protein [SPL), which 
are involved in leaf and flower development. The miRl 60 target 
mRNA was annotated with an auxin response factor {ARFj, which 
might be associated with auxin response during plant organogen- 
esis [28]. These results provide valuable information for the 
fiirther research. 

To describe the gene functions, we classified the potential 
targets into three categories based on TAIR GO annotations: 
molecular functions, biological processes and cellular components. 
For the molecular function category, genes were assigned to eight 
subcategories (Figure 5A). The GO terms binding (35.06%), 
enzyme regulator activity (22.87%), and transcription factor 
activity (14.33%) were the most frequent, and especially RNA 
polymerase II transcription cofactor activity, which belonged to 
the transcription factor activity subcategory, is important for 



miRNA regulation [29]. Biological processes, which included 11 
subcategories (Figure 5B), most frequently included GO terms 
involved in developmental processes (14.87%) and regulation of 
transcription (11.21%). Five subcategories of cellular components 
were identified (Figure 5C), of which the most frequent were 
nucleus (20.64%) and membrane (12.22%). Thus, the high 
frequency of GO terms associated with developmental processes 
indicated that many of the miRNAs identified in this study were 
involved in A. crassicarpa organogenesis by regulating molecular 
functions, biological processes and cellular components. 

Expression Patterns of miRNAs During Organogenesis 

To further investigate the role of miRNAs during A. crassicarpa 
organogenesis, 14 miRNAs of known function or high expression 
counts (one novel and 13 conserved [30-43]) (Table 6) were 
selected for analysis by qRT-PCR [44] , which is a reliable method 



Table 4. Novel mIRNA families identified from Acacia crassicarpa. 



mlRNA 


Location 


Sequence (5 -3 ) 


Length 


Count 


MFEI 


acr-novell 


gi|293578058:l:181- 


AUUAUAGGAACACUUUUGUAG 


21 


22 


1.6 


acr-novell* 




AAAAGUGUUCCUAUAAUUAGGACC 


24 


2 


1.6 


acr-novel2 


gi|293579709:l:247-l- 


UGUGGGACGCUAGGAAAGAG 


20 


2 


0.9 


acr-novel2* 




UCUUUCCUAUGCCUCCCAUACC 


22 


4515 


0.9 


acr-nove!3 


gi|293581 175:360:620- 


UCCGGUCUUAAUUAUAAGAAC 


21 


7 


1.6 


acr-novel4 


gma_Gm07:35780857:35781 1 1 8-1- 


AUCGACGCUGCACUCAAUCAUU 


22 


58 


1 


acr-novel5 


ptr_Chrl 3:84961 56:849641 6- 


UACUUCAAGUAGAGCAACCCC 


21 


7 


0.9 



doi:l 0.1 371 /journal.pone.0093438.t004 



PLOS ONE I www.plosone.org 10 April 2014 | Volume 9 | Issue 4 | e93438 



miRNAs during Organogenesis of Acacia crassicarpa 



Table 5. Predicted target genes of mIRNAs in Acada crassicarpa. 





miRNA 


Target 


Score 


Annotation 


MIR156 


AT1G69170 


1 


SBP domain transcription factor family protein 




AT5G50570 


1 


SBP domain transcription factor family protein 




AT5G50670 


1 


SBP domain transcription factor family protein 




AT1G53160 


2 


squamosa promoter binding proteln-IIke 4 




ATI G7 1690 


3 


Protein of unknown function (DUF579) 




AT3G15270 


3 


squamosa promoter binding protein-like 5 




AT5G 10945 


3 


MIR156D; miRNA 




AT5G 10946 


3 


unknown protein 




AT2G 19420 


0 


unknown protein 




ATI G29900 


2 


carbamoyl phosphate synthetase B 




AT2G33810 


2 


squamosa promoter binding proteln-IIke 3 




AT3G28690 


2 


Protein kinase superfamily protein 




AT1G17500 


3 


ATPase E1-E2 type family protein/haloacid 
dehalogenase-like hydrolase family protein 




AT1G29030 


3 


Apoptosis Inhibitory protein 5 (APIS) 




ATI G44790 


3 


ChaC-IIke family protein 




AT1G48410 


3 


Stabilizer of iron transporter SufD/Polynucleotldyl 
transferase 




AT1G61580 


3 


R-protein L3 B 




ATI G62305 


3 


Core-2/l-branching 

beta-1,5-N-acetylglucosaminyltransferase 
family protein 




ATI G7 1240 


3 


Plant protein of unknown function (DUF639) 




AT2G11910 


3 


unknown protein 




AT3G04870 


3 


zeta-carotene desaturase 




AT3G06670 


3 


binding 




AT3G18210 


3 


2-oxoglutarate {20G) and Fe(ll)-dependent 
oxygenase superfamily protein 




AT5G32566 


3 


transposable element gene 




AT1G27360 


0 


squamosa promoter-like 1 1 




AT1G27370 


0 


squamosa promoter binding proteln-IIke 10 




AT2G42200 


0 


squamosa promoter binding protein-like 9 




AT3G57920 


0 


squamosa promoter binding proteln-like 15 




AT5G43270 


0 


squamosa promoter binding protein-like 2 




ATI G22000 


3 


FBD, F-box and Leucine Rich Repeat 
domains containing protein 




AT3G44096 


3 


transposable element gene 




AT3G47170 


3 


HXXXD-type acyl-transferase family protein 




AT5G66600 


3 


Protein of unknown function, DUF547 




AT1G16916 


3 


unknown protein 




AT4G24270 


3 


EMBRYO DEFECTIVE 140 




AT5G11977 


3 


MIR156E; miRNA 




AT5G08620 


3 


DEA(D/H)-box RNA helicase family protein 




AT5G24930 


3 


CONSTANS-like 4 




AT5G44280 


3 


RING lA 




AT5G55590 


3 


Pectin lyase-IIke superfamily protein 




AT5G55835 


3 


MIR156H; miRNA 




AT5G20380 


3 


phosphate transporter 4;5 




ATI G48742 


2 


MIR157D; miRNA 




AT3G18217 


2 


MIR157C; miRNA 
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miRNA 


Target 
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ATI G30450 


3 


cation-chloride co-transporter 1 




ATI G73650 


3 


Protein of unknown function (DUF1295) 




AT3G 13490 


3 


Lysyl-tRNA synthetase, class II 




AT4G31050 


3 


Biotin/lipoate A/B protein ligase family 




AT5G35930 


3 


AMP-dependent synthetase and ligase 
family protein 




AT4G05400 


3 


copper ion binding 


MIR1 59 


AT2G3401 0 


3 


unknown protein 




AT4G37770 


2 


1-amino-cyclopropane-1-carboxylate synthase 8 




AT3G60460 


3 


myb-like HTH transcriptional regulator 
family protein 




AT4G2241 5 


3 


transposable element gene 




AT3G2891 5 


3 


transposable element gene 




AT3G33133 


3 


transposable element gene 




AT5G65540 


3 


unknown protein 




AT2G26950 


2 


myb domain protein 104 




AT3G1 1440 


2 


myb domain protein 65 




AT3G52030 


2 


F-box family protein with WD40/YVTN 
repeat doamin 




AT5G06100 


2 


myb domain protein 33 




ATI G08030 


3 


tyrosyl protein sulfotransf erase 




AT1G52100 


3 


Mannose-binding lectin superfamily protein 




ATI G76390 


3 


ARM repeat superfamily protein 




AT2G32460 


3 


myb domain protein 101 




AT3G33076 


3 


transposable element gene 




AT3G33084 


3 


transposable element gene 




AT4G26930 


3 


myb domain protein 97 




AT4G27330 


3 


sporocyteless (SPL) 




AT5G16810 


3 


Protein kinase superfamily protein 




AT5G28335 


3 


transposable element gene 




AT5G55020 


3 


myb domain protein 120 


MIR160 


AT2G28350 


2 


auxin response factor 10 


MIR164 


AT1G56010 


2 


NAC domain containing protein 1 




AT3G15170 


3 


NAC domain transcriptional regulator 
superfamily protein 




AT5G53950 


3 


NAC domain transcriptional regulator 
superfamily protein 


MIR167 


AT3G04765 


3 


MIR167C; miRNA 




AT3G07120 


3 


RING/U-box superfamily protein 


MIR169 


AT3G05690 


2 


nuclear factor Y, subunit A2 




AT5G06510 


2 


nuclear factor Y, subunit AID 




AT1G68560 


3 


alpha-xylosidase 1 




AT2G46020 


3 


transcription regulatory protein SNF2, 
putative 




A 1 DVjZd I 1 U 


3 


z.inc linger, t^^nt-^t Type (kiinij nngerj 
family protein 




ATI G 17590 


3 


nuclear factor Y, subunit A8 




AT1G54160 


3 


nuclear factor Y, subunit A5 




AT5G42120 


3 


Concanavalin A-like lectin protein kinase 
family protein 




AT5G 12840 


2 


nuclear factor Y, subunit A1 




AT3G20910 


3 


nuclear factor Y, subunit A9 
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MIR171 


AT2G45160 


1 


GRAS family transcription factor 




AT3G60630 


1 


GRAS family transcription factor 




AT4G00150 


1 


GRAS family transcription factor 




AT1G62035 


3 


MIR171C; miRNA 


MIR172 


AT5G60120 


1 


target of early activation tagged (EAT) 2 




ATS G6 5 790 


2 


myb domain protein 68 




AT5G67180 


2 


target of early activation tagged (EAT) 3 




ATI G6 1290 


3 


syntaxin of plants 124 




AT3G 11435 


3 


MIR172C; miRNA 




AT5G 12900 


3 


unl(nown protein 




AT4G36920 


1 


Integrase-type DNA-binding superfamily 
protein 




AT2G28056 


2 


MIR172/MIR172A; miRNA 




AT2G28550 


2 


related to AP2.7 




AT2G39250 


3 


Integrase-type DNA-binding superfamily 
protein 




AT3G49690 


3 


myb domain protein 84 




AT3G 54990 


2 


Integrase-type DNA-binding superfamily 
protein 




AT3G55512 


3 


MIKI 72D; miRNA 




AT5G04275 


3 


MIRI 72/(vilRl 72b; miRNA 




AT5G59505 


3 


MIRI 72E; miRNA 


MIR31 9 


AT2G26950 


2 


myb domain protein 104 




AT3G1 1440 


2 


myb domain protein 65 




AT5G06100 


2 


myb domain protein 33 




AT1G53230 


3 


TEOSINTE BRANCHED 1, cycloidea and 
PCF transcription factor 3 




AT3G33076 


3 


transposable element gene 




AT3G33084 


3 


transposable element gene 




AT3G66658 


3 


aldehyde dehydrogenase 22A1 




AT4G2241 5 


3 


transposable element gene 


MIR390 


AT2G38325 


3 


MIR390A; miRNA 


MIR393 


AT3G62980 


1 


F-box/RNI-like superfamily protein 




AT4G03 1 90 


2 


GRRI-like protein 1 


MIR394 


ATI G27340 


1 


Galactose oxldase/kelch repeat superfamily 
protein 




AT3G29660 


3 


transposable element gene 




AT5G09670 


3 


loricrin-related 




AT5G09672 


3 


conserved peptide upstream open reading frame 21 


MIR395 


ATI G26973 


2 


MIR395A; miRNA 




AT1G69792 


2 


MIR395D; miRNA 




ATI G69795 


2 


MIR395E; miRNA 




AT1G26975 


3 


MIR395B; miRNA 




AT1G26985 


3 


MIR395C; miRNA 




ATI G69797 


3 


MIR395F; miRNA 




AT2G41060 


3 


RNA-binding (RRM/RBD/RNP motifs) 
family protein 




AT3G32280 


3 


ATP-dependent helicase family protein 




AT4G 14770 


3 


TESMIN/TSOI-like CXC 2 




AT5G04140 


3 


glutamate synthase 1 
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AT5G 19900 


3 


PRLI-interacting factor, putative 




AT2G28780 


2 


unknown protein 




ATI G50930 


1 


unknown protein 




A 1 Dtj4J/oU 


1 


Pseudourldine synthase/archaeosine 
transglycosylase-like family protein 




AT3G22890 


3 


ATP sulfurylase 1 




AT4G 14680 


3 


Pseudouridine synthase/archaeosine 
transglycosylase-like family protein 




AT4G06509 


3 


transposable element gene 




AT5G 13630 


3 


magnesium-chelatase subunit chlH, 
chloroplast, putative/Mg-protoporphyrin 
IX chelatase 


MIR396 


AT5G14130 


3 


Peroxidase superfamily protein 




AT1G58020 


2 


transposable element gene 




AT2G45480 


3 


growth-regulating factor 9 




AT3G33106 


3 


transposable element gene 




AT5G01370 


3 


ALC-interacting protein 1 




AT5G07700 


3 


myb domain protein 76 




AT5G35407 


1 


MIR396B; miRNA 




AT2G36400 


3 


growth-regulating factor 3 




AT3G01910 


3 


sulfite oxidase 




AT2G04270 


3 


RNAse E/G-like 




AT3G52910 


3 


growth-regulating factor 4 




AT5G 16690 


3 


origin recognition complex subunit 3 


MIR397 


AT5G 18420 


3 


unknown protein 




AT2G29130 


2 


laccase 2 




A 1 4LiUd I Ud 


2 


IVIIKjy/A, minlNA 


jviiKjyes 


A 1 DK3 \ 43dU 


1 


Core-2/l-branching 

beta-1 ,5-N-acetylglucosaminyl 

transferase family protein 


MIR403 


ATI G3 1280 


0 


Argonaute family protein 




ATI G03060 


3 


Beige/BEACH domain;WD domain, 
G-beta repeat protein 




AT4G30620 


3 


Uncharacterised BCR, YbaB family 
COG0718 


MIR408 


AT2G47020 


0 


Peptide chain release factor 1 




AT2G02850 


1 


plantacyanin 


MIR857 


ATI G 14610 


2 


valyl-tRNA synthetase/valine-tRNA 
ligase (VALRS) 




ATI G 74060 


2 


Ribosomal protein L6 family protein 




AT5G28350 


2 


Quinoprotein amine dehydrogenase 




ATI G06670 


3 


nuclear DEIH-boxhelicase 




ATI G 17040 


3 


SH2 domain protein A 




ATI G33475 


3 


SNARE-like superfamily protein 




ATI G43590 


3 


transposable element gene 




AT1G48010 


3 


Plant invertase/pectin methylesterase 
inhibitor superfamily protein 




AT3G20090 


3 


cytochrome P450, family 705, subfamily A, 
polypeptide 18 




AT3G54230 


3 


suppressor of abi3-5 




AT3G57660 


3 


nuclear RNA polymerase A1 



PLOS ONE I www.plosone.org 14 April 2014 | Volume 9 | Issue 4 | e93438 



miRNAs during Organogenesis of Acacia crassicarpa 



Table 5. Cont. 





miRNA 


Target 


Score 


Annotation 




AT3G61480 


3 


Quinoprotein amine dehydrogenase, beta 
chain-like; RICI-like guanyl-nucleotide 
exchange factor 




AT4G08460 


3 


Protein of unknown function (DUF1644) 




AT4G10170 


3 


SNARE-like superfamily protein 




AT5G06040 


3 


self-incompatibility protein-related 




AT5G59470 


3 


Mannose-P-dolichol utilization defect 1 protein 


MIR858 


AT1G06180 


3 


myb domain protein 13 




ATI G66230 


3 


myb domain protein 20 




AT2G47460 


3 


myb domain protein 12 




AT3G08500 


3 


myb domain protein 83 




AT4G12350 


3 


myb domain protein 42 




AT5G35550 


3 


Duplicated homeodomain-like 
superfamily protein 




AT5G49330 


3 


myb domain protein 111 


MIR1310 


AT2G 16592 


3 


Bifunctional inhibitor/lipid-transfer protein/seed 
storage 2S albumin superfamily protein 




AT3G29050 


3 


receptor-like protein kinase-related 


MIR1448 


AT5G40100 


1 


Disease resistance protein (TIR-NBS-LRR class) 
family 




AT1G65790 


2 


receptor kinase 1 




ATI G65800 


2 


receptor kinase 2 




AT3G46530 


2 


NB-ARC domain-containing disease 
resistance protein 




AT1G63350 


3 


Disease resistance protein {CC-NBS-LRR class) 
family 




AT2G30740 


3 


Protein kinase superfamily protein 




AT4G00420 


3 


Double-stranded RNA-binding domain 
(DsRBD)-containing protein 


MIR1 51 1 


AT3G28160 


2 


transposable element gene 




AT5G42820 


2 


Zinc finger C-x8-C-x5-C-x3-H type family protein 




ATI G03050 


3 


ENTH/ANTH/vl-lb superfamily protein 




ATI G 17340 


3 


Phosphoinositide phosphatase family protein 




AT1G61610 


3 


S-locus lectin protein kinase family protein 




AT3G21360 


3 


2-oxoglutarate (20G) and Fe(ll)-dependent 
oxygenase superfamily protein 




AT4G24450 


3 


phosphoglucan, water dikinase 




AT5G27965 


3 


transposable element gene 




AT5G32690 


3 


Pseudogene of AT2G29880 


MIR2111 


AT3G27150 


3 


Galactose oxidase/kelch repeat superfamily protein 


MIR2592 


AT4G11040 


3 


Protein phosphatase 2C family protein 


MIR2916 


ATI G 15780 


3 


unknown protein 




AT1G50310 


3 


sugar transporter 9 




AT1G65875 


3 


pseudogene 




AT1G66120 


3 


AMP-dependent synthetase and ligase family protein 






3 


RNA polymerase II transcription mediators 




AT5G66710 


3 


Protein kinase superfamily protein 


MIR4414 


AT2G04860 


3 


Tetratricopeptide repeat (TPR)-like 
superfamily protein 


MIR4415 


AT3G 14490 


3 


Terpenoid cyclases/ProteIn prenyltransferases 
superfamily protein 


MIR5139 


AT3G28160 


0 


transposable element gene 




AT5G27965 


1 


transposable element gene 
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AT4G05360 


2 


Zinc knuckle (CCHC-type) family protein 




AT5G67460 


2 


0-Glycosyl hydrolases family 17 protein 




ATI G34470 


3 


Protein of unknown function (DUF803) 




ATI G64470 


3 


Ubiquitin-like superfamlly protein 




ATI G72800 


3 


RNA-binding (RRM/RBD/RNP motifs) 
family protein 




AT2G 26430 


3 


arginine-rich cyclin 1 




AT3G 50650 


3 


GRAS family transcription factor 




AT3G53270 


3 


Small nuclear RNA activating complex 
{SNAPc}, subunit SNAP43 protein 




AT4G09880 


3 


unknown protein 




AT5G34450 


3 


transposable element gene 


MIR5256 


AT3G47990 


3 


SUGAR-INSENSITIVE 3 




AT5G47420 


3 


Tryptophan RNA-binding attenuator protein-like 




AT5G 54690 


3 


galacturonosyltransferase 12 


MIR5368 


AT5G5 1 1 30 


3 


S-adenosyl-L-methionine-dependent 
methyltransferases superfamlly protein 




AT4G33920 


3 


Protein phosphatase 2C family protein 


MIR5740 


AT5G44925 


3 


transposable element gene 




AT5G50180 


3 


Protein kinase superfamlly protein 


MIR6478 


AT4G22270 


0 


Protein of unknown function (DUF3537) 


acr-novel2* 


AT4G 19520 


1 


disease resistance protein {TIR-NBS-LRR class} family 




AT4G09430 


2 


Disease resistance protein {TIR-NBS-LRR class) family 




AT4G37190 


3 


LOCATED IN: cytosol, plasma membrane; 




AT5G45260 


3 


Disease resistance protein (TIR-NBS-LRR class) 



doi:1 0.1 371 /journal.pone.0093438.t005 



to detect and measure the expression levels of miRNAs in various 
tissues. The relative quantitative results demonstrated that all 
selected miRNAs were expressed in different organogenetic tissues 
and with entirely different expression profiles. On the basis of 
abundance trends at the six developmental stages, the 14 miRNAs 
were divided into four clusters through multivariate statistical 
analysis using general cluster analysis procedure of SPSS Statistics 
18.0 (Figure 6). 

The first cluster, comprising acr-miR159a, acr~miR3 1 9a, acr- 
miR162, acr-miRl 7 1 q, acr-miR390a, and acr-miR396, exhibited 
similar expression patterns and showed a low expression level 
compared with the other three groups. These miRNAs were 
barely accumulated during stages S 1 and S2 and were induced at 
stages S3 and S4 when the peak expression level was observed 
(Figure 7A-F), which implied that their accumulation may be 
required for embryogenic callus formation. The second cluster, 
consisting of acr-miRl 64a, acr-miRl 67a, and acr-miRl 68a, showed a 
stage-specific expression pattern. The expression level peaked 
during S4 but at the other stages expression was less abundant or 
undetectable (Figure 7G-I). The cultures were transferred to light 
at the S4 stage, which suggested these miRNAs might have an 
important relationship with light. In the third cluster, acr-miR156, 
acr-miRl 66, and acr-novel2* were expressed at varied levels in the 
different stages. Their major peak in expression was observed at 
S3, whereas expression at the other stages was relatively lower 
(Figure 7J-L). The striking differences in expression inferred their 
fijnction is exerted in different developmental stages. The 
remaining two miRNAs {acr-miR397 and acr-miR39S) showed 



relatively high expression levels in most of the detected tissues, 
especially acr-miR397. In this group, the miRNAs expression level 
rose continuously during the successive developmental stages 
except at S4, and accumulated in S3 and reached their highest 
levels in the adventitious shoots at S6 (Figure 7M-N), which 
suggested these miRNAs play a major regulatory role in 
morphogenesis during advanced differentiation of A. crassicarpa 
adventitious shoots. 

The above results concluded that the different miRNAs were 
indicated to have functions at different developmental stages 
during A. crassicarpa organogenesis. 

Discussion 

In vitro organogenesis refers to the process involving regenera- 
tion of adventitious shoots, adventitious roots or other organs from 
plant tissue or cell aggregates (callus), which is an important aspect 
of plant development [45] . Many genes involved in the regulation 
of organogenesis have been investigated [46], but the molecular 
mechanisms underlying the process are stiU not well understood. 
Previous studies have confirmed that miRNAs play important 
roles in a variety of developmental processes. However, the role of 
miRNAs during organogenesis is poorly studied. In the present 
investigation, we used a high-throughput strategy to perform large- 
scale cloning and characterization of miRNAs involved in A. 
crassicarpa adventitious shoot organogenesis. We identified 189 
known and 7 novel miRNAs from more than 57 miRNA families, 
and 1 4 conserved and novel miRNAs were selected for analysis of 
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Figure 5. Functional classification of microRNA targets. Gene ontology of the predicted targets for 37 differentially expressed miRNAs. (A) 
Molecular functions was divided into 10 functional groups. (B) Biological processes was divided into 14 functional groups. (C) Cellular components 
was divided into seven functional groups. 
doi:1 0.1 371 /journal.pone.0093438.g005 



expression patterns by qRT-PCR. These results provide valuable 
information on the molecular mechanism of organogenesis in A. 
crassicarpa. Furthermore, target gene prediction and GO annota- 
tion demonstrated that putative miRNA targets were involved in a 
broad variety of regulatory events, including molecular functions, 
biological processes and cellular components. According to specific 
stages, the expression patterns of 14 selected miRNAs were 
observed. 

Micropropagation is an important and reliable technique for the 
production of large quantities of many plant species, particularly as 
a tool for large-scale plant breeding programs, and is used as a 
model system to research cell differentiation and development 
during organogenesis [47] . Changes in miRNA expression during 
plant somatic embryogenesis has been confirmed in many plants 
species, such as A. thaliana [48], rice [49], Japanese larch [50], 
Valencia sweet orange [51], Liriodendron chinense [52], loblolly pine 
[53], and longan [54]. However, the function of miRNAs during 
plant organogenesis is poorly studied. In order to better 
understand the potential modulation of miRNAs in organogenesis, 
we studied the expression patterns of miRNAs utilizing the A. 
crassicarpa organogenesis in vitro culture system described by Yao 
et al [18]. Combined with a bioinformatic analysis and comparison 
with miRNA families that are expressed during somatic embryo- 
genesis in other plant species, the results revealed some of the 
important miRNAs functions related to plant regeneration. 
Dilferences in expression pattern were observed among 14 
miRNAs analyzed by qRT-PCR. One group of miRNAs [Acr- 
miR159a, acr-miR3 1 9a, acr-miR162, and acr-miR171(i) were up- 
regulated at S3 and the expression level peaked at S4 (Figure 7). 
MiR159 is a widespread gene present in many plant and animal 
species [55] . As a highly conserved miRNA in the plant kingdom 
miR159 negatively controls gene expression by targeting mainly 



MYB33 and MYB65 during seed germination and floral 
development [56-57]. In the present study, acr-miR159a transcripts 
accumulated until embryogenic callus was present at stage S3 of ^. 
crassicarpa organogenesis, which suggested that acr-miR159a might 
modulate gene expression during embryogenic callus differentia- 
tion. This result is similar to the pattern observed in larch [50] and 
longan [54]. Acr-miR319a, which has an extremely similar 
sequence to that of the miR159 family, showed similar relative 
expression levels during A. crassicarpa organogenesis [58]. Acr- 
miR319a and acr-miR159, showed the same putative target (MYB 
transcription factors) in the regulation of organogenesis in die 
present study, which suggests they play similar regulatory roles. 
The functions of miR171 and miR162 are rarely reported. In A. 
thaliana, the potential target of miRl 71 is the scarecrow-like [SCE) 
gene family possessing the GRAS domain [59], which is consistent 
with our findings (Table 5). In our study, these two miRNAs were 
up-regulated at S3 and S4, indicating that miRl 71 and miR162 
might function in the embryogenic callus, which is consistent with 
rice somatic embryogenesis [49]. Expression of the remaining 
miRNAs of the first cluster, acr-miR390a and acr-miR396, was 
elevated at S3 then declined, but was maintained at a relatively 
high level from S3 to S5. MiR396 is reported to target growth- 
regulating factor [GRF] transcription factors, performing a negative 
coordinating role in leaf cell proliferation of A. thaliana [60]. 
Overall, members of the first cluster might play regulatory roles in 
the embryogenic callus stages oi A. crassicarpa organogenesis. 

S4 is a crucial stage during A. crassicarpa organogenesis when 
buds are induced and began to turn green, and are ready to 
develop into adventitious shoots. miRNAs of the second cluster, 
consisting of acr-miRl 64a, acr-miRl 67a, and acr-miRl 68a, exhibited 
stage-specific expression, which indicated that they may have 
stage-specific functions during bud formation. These three 



Table 6. miRNAs associated with plant development. 
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Figure 6. Cluster analysis of expression trends of selected miRNAs. Dendrogram based on expression trends of 14 selected miRNAs during 
the entire Acacia crassicarpa organogenesis process. 
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miRNAs were expressed at low or undetectable levels at aJl stages 
except S4, during which expression peaked. In Arabidopsis, miRl 64 
mainly controls the NAM/ATAF/CUC {MACl) domain-tran- 
scription factor family. NACl is involved in transitions in auxin 
signalling, and facilitates growth of lateral roots [61]. CUC 
regulates organ separation from fasciculate buds during embryo- 
genesis. In the present study, the same putative target (NAC- 
domain transcription factor) implied that miRl 64 has the same 
function in A. crassicarpa. MiRl 68 targets the AGOl gene, which is 
involved in plant development by feed-back regulation. The 
expression patterns of acr-miRl 64a and acr-miRl 68a were consis- 
tent with results reported for citrus after callus is cultured in the 
light [51]. However, our findings were inconsistent with the 
pattern observed in larch, in which miRl 68 was not expressed at a 
notably high level when callus began to turn green [50]. These 
conflicting results suggest that the function of miRl 68 might vary 
between species. Aux/IAA and ART, the putative targets of 
miRl 67, are two important protein families that respond to auxin 
signal. Ru et al. [62] suggested that the high abundance of ARF8 
led to the low expression level of miRl 67. These results imply the 
target of miRl 67 conservative property. In conclusion, this cluster 
of miRNAs may function in the bud formation stages and response 
to light in A. crassicarpa organogenesis. 

The important events that occur during A. crassicarpa adventi- 
tious shoot organogenesis from a zygotic embryo involve 
redifferentiation. In rice [49], miRl 56 is important in the transition 
from undifferentiated to differentiated callus during somatic 
embryogenesis by targeting SPL genes. Zhang et al. [50] and 
Wu et al. [51] reported that SPL controls the somatic embryo 
induction process. SPL genes also have other regulatory roles in 
different biological processes, such as induction of the floral 
transition and consequent shortening of the vegetative phase [63], 
and regulation of the juvenile to adult transition during plant 
development [64]. In the present study, acr-miRl 5 6r was accumu- 



lated at S3 while the zygotic embryo is differentiating into 
embryogenic callus during organogenesis. As shown in Figure 7, 
the abundance of the other members of the third cluster {acr- 
miRl 66h and acr-novel2*) also increased continuously until S3,but 
exhibited very low or undetectable expression levels during S4— S6 
and especially at S5. Given that these three miRNAs showed 
similar expression patterns, we hypothesize that they modulate 
redifferentiation during induction of organogenesis in A. crassicarpa 
as the fatal genes. 

With regard to the fourth cluster, relatively high expression of 
acr-miR397 and acr-miR39d was observed. MiR397, for which a 
target gene of laccase is implicated, is associated with lignin 
biosynthesis and primary cell wall [65-66]. The detection of 
miR397 and corresponding targets occurs during the regulation of 
lignification and thickening of the cell wall in secondary cell 
growth in rice, larch and citrus [49-51]. MiR398 is known to 
target Cu/Zn superoxide dismutases (CDSs), which are associated 
with stress response [43] . In the present study, the levels of acr- 
miR397 and acr-miR398 increase continually until S6, at which 
stage the expression level peaked, except for a slight decline at S4. 
Zygotic embryos undergo differentiation at stages SI -S3 until the 
formation of buds at S4, which suggests that laccase, the target 
gene of miR397, regulates lignification and cell wall thickness 
during organogenesis. The decrease in expression level of miR397 
at S4, when a cluster of adventitious shoots had developed, is 
attributed to negative modulation of the formation of thickened 
cell walls for adventitious shoots. 

In summary, a global analysis of miRNAs expression during A. 
crassicarpa adventitious shoot organogenesis was carried out. The 
results of a bioinformatic analysis and experimental tests revealed 
putative regulatory functions for the miRNAs in Acacia crassicarpa 
organogenesis. These findings provide important information for 
deep sequencing research of miRNAs and future large-scale 
propagation and breeding of leguminous trees. 
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Figure 7. qRT-PCR analysis of relative expression levels of selected miRNAs at six stages of Acacia crassicarpa organogenesis. The fold 
change in gene expression was normalized to controls (mature zygotic embryo) with the 2~ method using 5.8S rRNA as an internal standard. 
Templates for all miRNAs real-time PCR were 1/20 dilutions of original cDNAs reverse-transcribed from 300 ng miRNA. Each bar shows the mean of 
triplicate assays. 
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